Many-body corrections to cyclotron resonance in monolayer and bilayer graphene 



O 

(N 

C 

a 

(N 



T3 

c 

o 

> 



o 



X 



K. Shizuya 
Yukawa Institute for Theoretical Physics 
Kyoto University, Kyoto 606-8502, Japan 

Cyclotron resonance in graphene is studied with focus on many-body corrections to the resonance 
energies, which evade Kohn's theorem. The genuine many-body corrections turn out to derive from 
vacuum polarization, specific to graphene, which diverges at short wavelengths. Special emphasis 
is placed on the need for renormalization, which allows one to determine many-body corrections 
uniquely from one resonance to another. For bilayer graphene, in particular, both intralayer and 
interlayer coupling strengths undergo infinite renormalization; as a result, the renormalized velocity 
and interlayer coupling strength run with the magnetic field. A comparison of theory with the 
experimental data is made for both monolayer and bilayer graphene. 

PACS numbers: 73.422.Pr,73.43.Lp,76.40.+b 



I. INTRODUCTION 

Graphene, a monolayer graphite, attracts great atten- 
tion for its unusual electronic transporter— as well as 
its potential applications. It supports as charge carri- 
ers massless Dirac fermions, which lead to such exotic 
phenomena as the half-integer quantum Hall effect and 
minimal conductivity. 

The multi-spinor character of the electrons in graphene 
derives from the sublattice structure of the underlying 
honeycomb lattice, and this immediately implies, in the 
low-energy effective theory of graphene, the quantum na- 
ture of the vacuum stated the conduction and valence 
bands are related by charge conjugation and, in partic- 
ular, the latter acts as the Dirac sea. Graphene in a 
magnetic field B thus gives rise to a particle-hole sym- 
metric "relativistic" pattern of Landau levels, with spec- 
tra e„ oc ± y/\n\\/B unequally spaced, together with four 
characteristic zero-energy Landau levels (whose presence 
has a topological origin 8 -). 

This nontrivial vacuum structure is the key feature 
that distinguishes graphene and its multilayers from con- 
ventional quantum Hall (QH) systems. In particular, bi- 
layer graphene 9 supports, as a result of interlayer cou- 
pling, massive fermions, which, in a magnetic field, again 
develop a particle-hole symmetric tower of Landau levels, 
with an octet of zero-energy levels^ Bilayer graphene 
has a unique property that its band gap is externally 
controllable*^ - — 

Graphene and its multilayers give rise to rich spec- 
tra of cyclotron resonance, with resonance energies vary- 
ing from one transition to another within the electron 
band or the hole band, and, notably, even between the 
two bands. This is in sharp contrast with conventional 
QH systems with a parabolic dispersion, where cyclotron 
resonance (optically-induced at zero momentum transfer 
k = 0) takes place between adjacent Landau levels, hence 
at a single frequency w c = eB/m* which, according to 
Kohn's theorem,— is unaffected by Coulomb interactions. 
Nonparabolicityi^ of the electronic spectra in graphene 
evades Kohn's theorem and offers the possibility to de- 



tect the many-body corrections to cyclotron resonance, 
as discussed theoretically for monolayer graphene ! 17 ' 18 

Experiment has already studied via infrared spec- 
troscopy cyclotron resonance in monolaye r 19 ' 20 and bi- 
layer— graphene, and verified the characteristic features 
of the associated Landau levels. Data generally show no 
clear sign of the many-body effect, except for onei 9 - on 
monolayer graphene. 

The purpose of this paper is to study the many-body 
effect on cyclotron resonance in graphene, by construct- 
ing an effective theory of cyclotron resonance within the 
single-mode approximation. It is shown that the gen- 
uine many-body corrections arise from vacuum polariza- 
tion, specific to graphene, which actually diverges loga- 
rithmically at short wavelengths and requires renormal- 
ization. Our approach in part recovers results of ear- 
lier studie s 17 ' 18 on monolayer graphene but essentially 
differs from them in this handling of cutoff-dependent 
corrections by renormalization, which allows one to de- 
termine many-body corrections uniquely from one res- 
onance to another. Our analysis also reveals that for 
bilayer graphene both intralayer and interlayer coupling 
strengths undergo renormalization! 22 ' 23 We compare the- 
ory with the experimental data for monolayer and bilayer 
graphene. 

In Sec. II we briefly review the effective theory of 
graphene and, in Sec. Ill, study cyclotron resonance in 
monolayer graphene, with focus on the many-body cor- 
rections and renormalization. In Sec. IV we extend our 
analysis to bilayer graphene. Section V is devoted to the 
summary and discussion. 



II. MONOLAYER GRAPHENE 

Graphene has a honeycomb lattice which consists of 
two triangle sublattices of carbon atoms. The electrons 
in graphene are described by a two-component spinor 
field (ipA,ipB) on two inequivalent lattice sites (A,B). 
The electrons acquire a linear spectrum near the two in- 
equivalent Fermi points (K and K') in the Brillouin zone, 
with the "light velocity" v = (V3/2) a L j /h w 10 6 m/s 
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related to the intralayer coupling 70 = jab ~ 2.9 eV 
(with a L = 0.246nm). 

Their low-energy features are described by an effective 
Hamiltonian of the form 7 



d 2 x 



U± — vo (crilli + cr 2 n 2 =F 5m(Tz) — eA Q , (2.1) 

where 11^ = —idi + eAi [i — (1,2) or (x,y)] involve 
coupling to external electromagnetic potentials A a = 
(Ai,A ). Here tp = (^ljV^)* = (^a^bY stands for the 
electron field at one (or K) valley, and ip' = (V'biV'a)* 
to one at another valley, with A and B referring to the 
associated lattice sites. 8m denotes a possible tiny asym- 
metry in sublattices. 

Let us place graphene in a uniform magnetic field B± = 
B > 0; we set Ai — > £>(— y,0). The electron spectrum 
then forms an infinite tower of Landau levels L n of energy 



= s n uj c y/\n\ + (Sm) 2 e 2 /2, 



(2.2) 



labeled by integers n — 0,±1,±2,..., and p x (or y = 
£ 2 p x with the magnetic length I = \j\feB). Here s n = 
sgn{n} = ±1 specifies the sign of the energy e„, and 



uj c = V2v /£ w 36.3 x i; [10 6 m/s] y/B\T} mcV. (2.3) 

is the basic cyclotron frequency; -u [10 6m / s ] stands for v a 
in units of 10 6 m/s and B[T] for a magnetic field in tesla. 

Suppose, without loss of generality, that 5m > 0. Then 
the n = level at the K valley has positive energy 
£o+ = vo5m > while the n = level at the K' valley has 
negative energy eo_ — —v$5m. In general, the spectra at 
the two valleys are related as e„|x = —e- n \K'- With the 
electron spin taken into account (and Zeeman splitting 
ignored for simplicity), each Landau level is thus fourfold 
degenerate, except for the doubly-degenerate n = 0± lev- 
els split in valley. With this feature in mind, we shall set 
the asymmetry 5m — > in what follows. 

The Coulomb interaction is written as 



H 



Coul 



p-pPp ■, 



(2.4) 



where p p is the Fourier transform of the electron density 
p — ip^ip + ip''tp' (here tp^tp, e -S-j is summed over spinor 
and spin indices); v p = 27ra/(eb|p|) is the Coulomb po- 
tential with the fine-structure constant a = e 2 /(47reo) ~ 
1/137 and the substrate dielectric constant eb. 

The Landau-level structure is made explicit by pass- 
ing to the \n,yo) basis, with the expansion ip(x,t) — 
J2n, yo ( x \ n >yo) i>n{ya,t)- (For conciseness, we shall only 
display the ip sector from now on.) The Hamiltonian H 
is thereby rewritten as 



H a = dyo^Tpliyoit) e n ip n (y ,t), 



(2.5) 



and the charge density p- p (t) = f c? 2 xe ip ' x ip^ip as^i 
00 00 

_ \ n kn \ „kn jykn 

P-p - P-p - 9 p k p ' 

k.n— — 00 k.n— — oc 

K n = Tp / dy ^t(yo,t)e^ r Myo,t), (2.6) 



where 7 P = e ^p 2 /4. r _ (r x ,r y ) = (ie 2 d/dy ,y ) stands 
for the center coordinate with uncertainty [7*3., r v ] = i£ 2 . 
The charge operators R p obey the algebra 25 

l H k > K p J= d e H k+ p -5 e K k+ P > 

(2.7) 

where k'p = k p — i kxp with kxp = k x p y — k y p x . This 
actually consists of two algebras associated with in- 
tralevel center-motion 2 ^ and interlevel mixing 2 ^ of elec- 
trons. 

The coefficient matrix is given by 

9p l = \b k b n {f p k ^ M ^ + s k s n f p k ^), (2.8) 
where b n — 1 for n ^ and bo = \/2; 

for k > n > 0, and f p k = (/^p) 1 "; p = p x -ip y - Note 
that (/p" are essentially the same at the two valleys, i.e., 

<?p n k' = 9p n \K and gf-\ K , = gf + \ K for |fc| > 1 and 
\n\ > 1; one simply needs to specify n = 0± accordingly. 



III. CYCLOTRON RESONANCE 

In this section we study cyclotron resonance in mono- 
layer graphene. Let us first note that the charge operator 
P-p = 9p n Rp l m Eq. (12.61) annihilates an electron at the 
nth level L n and creates one at the kth level L k ■ One may 
thus associate it with the interlevel transition L n — > L k 
or regard it as an interpolating operator for the exciton 
consisting of a hole at L n and an electron at L k ■ 

To describe such inter-Landau-level excitations one can 
make use of a nonlinear realization of the algebra, 
as is familiar from the theory of quantum Hall ferromag- 
neti^l One may start with a given ground state \G) and 
describe an excited collective state \G) over it as a local 
rotation in (n, yo) space, 



\G) 



-iO 



\G), 



where the operator e with 



-1 ^fcn n, 

Ip ^p *>-p 



kn 



(3.1) 



(3.2) 



locally rotates \G) by "angles" $p™, which define textures 
in (n, yo) space. (Remember here that p p ~ R p n are 



3 



diagonal in spin so that carry the spin index as well, 
though it is suppressed. In principle, one has to retain in 
O all possible pairs and = 3>™p contributing 

to the L n — > Lk transition.) 

Repeated use of the charge algebra (|2 . T[) allows 
one to express the texture-state energy (G\H\G) — 
(G\e l °He- w \G) with H = H a + H Goul as a functional of 
<j>fcn Qr x _ S p ace representative <fr kn (x, t) . The kinetic 
term for <£> fcn is supplied from the electron kinetic term, 
and one can write the effective Lagrangian for $ as 

= (G\(id t -H)\G) = (G\e l ° {id t -H)e' l °\G). (3.3) 

This representation systematizes the single-mode approx- 
imation 25 (SMA) within a variational framework^ The 
present theory thus embodies the nonperturbative fea- 
tures of the SMA. 

Indeed, for a transition from the filled level L n (with 
density p n ) to the empty level Lf., one finds 



£• = Pn £ $ -p( i5 * - £ p XC ) $ p" + ■ ' ' • 



(3.4) 



with the excitation spectrum given by the SMA formula, 
4" = ( G H k > [H,p%]]\G)/{G\pl k p%\G), (3.5) 

i.e., as the oscillator strength divided by the static struc- 
ture factor, 



(G\p^p%\G}/n = p n ^\gZ 



kn\2 



(3.6) 



where ft = j e? 2 x. 

Let us first consider the n = 0± — > n = 1 transitions 
at filling factor v = 2, with all n < 0± levels filled and 
all n > 1 levels empty. A laborious direct calculation of 
ei - e + Ae^ with 



Eq. J33J yields e° xc 



A4° = Pov k ti\glV 9 k y+Y,wll P ,u, 

p 

' P ,k = EW'-is-pi 2 )-^^™. ( 3 - 7 ) 



n<0 



where 7^ = e~ e k I" 1 and c p k = cos(£ 2 pxk). 

When the n — 0± level is partially filled, Ae^° involves 
the following contribution, 



^1° = £^k7#l 2 s(P + k)/7k k 



+ (c P ,k3 p 1 /.^ -l)S(p)/7^ 



(3.8) 



where s(p) stands for the projected static structure fac- 
tor defined as (G\p°° Pp °\G)/n = p {p S P;0 + s(p)}. The 
determination of s(p) is a highly nontrivial task which re- 
quires an exact diagonalization study of model systems^ 
and is not attempted here. We instead focus on the case 
of integer filling, for which s(p) is taken to vanish. 

Equation (I3.7[) . together with Eq. (I3.8p . essentially 
agrees with the result of MacDonald and Zhang 2 ^ for 



the Lo — >• L\ transition in the standard QH system. The 
key difference is that quantum fluctuations in I p k now 
involve a sum Yl n <-i(' ' ') ovcr infinitely many Landau 
levels in the valence band (or the Dirac sea). Actually, 
one can verify that the SMA expressions (|3.7[) and (|3.8[) 
equally apply to a general interlevel transition L a —> Lb if 
one sets the superscripts — > a and 1 —¥ b, in an obvious 
fashion, and takes the sum J2 n over filled levels. 

To be precise, the — > 1 transition of our interest 
consists of four channels, (0+ — > 1)\k and (0_ — > 1)\k' 
each with spin s z — ±1/2. One therefore has to consider 
mixing of four $ 10 to determine e k xc . Actually only the 
first term oc «k l5k°l 2 m ^l- (E2J), which comes from the 
direct Coulomb exchange, is responsible for such mix- 
in g 17 i 18 because the rest of terms are diagonal in spin 
and valley. Such direct terms cx «k Iff'l'kl 2 (with b 7^ a) 
in general vanish for k — > 0, and mixing thus takes place 
only for k ^ 0. In what follows we focus on the k = 
excitation energies e^O' with no mixing taken into ac- 
count. In addition, we make no distinction between the 
0± levels because g p n are essentially the same at the two 
valleys, as noted in Sec. II. 

The cyclotron-resonance energy for a general L a — > Li, 
transition with the Landau levels filled up to n — rif is 
written as eg? Q = e& — e a 



Ae£to with 



Ae b 



k=0 



p 



\g-l\ 2 



9*9-1 



n<rif 



(3.9) 

This is the basic formula we use in what follows. Note 
that nf = 1, 0+, — 1, — 2, — 3, ... correspond to the filling- 
factors v — 4 nf + 2 = 6, 2, —2, —6, — 10, respectively. 
The X)n<n f (ISp n | 2 — Iff-^! 2 ) term refers to the change in 
quantum fluctuations, via the a — > b transition, of the 
filled states. Its structure is easy to interpret physically: 
As an electron is excited from L a to Lb, \n — a, yo) —¥ 
\n = b,y' ), virtual transitions from any filled levels to 
the \b, y'o) state are forbidden while those to the newly 
unoccupied \a, yo) state are allowed to start. 

For standard QH systems this correction Ae^J van- 
ishes for each transition to the adjacent level, L n — > 
according to Kohn's theorem^ Indeed, one can 
verify, for the — > 1 transition, the relation 



\9 



00 P | 2 -I.9p 



10 1 2 



(3.10) 

£p/y2) and analogous ones 



(with g p ° -> 1 and g p ° - 
for other L n — > L n+1 as well. 

Interestingly, it happens that Eq. p.lOj) also holds 
for the 0± — > 1 transition in graphene, with g p = 1, 
g p ° = —£p/2 and gp 1 = 1 — £ 2 p 2 /4. Any nonzero shift 
Ae k =o f° r the — > 1 cyclotron resonance therefore comes 
from the quantum fluctuations of the Dirac sea, and ac- 
tually diverges logarithmically with the number Al of 
filled Landau levels in the sea, 



a4-° = E«p7 2 E M 

p -JV L <n<-l 

C N w (V2/8)(logA^- 1.017), 



( " 12 -\gi n \ 2 ) = v c c N , 



(3.11) 
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where 



V c = a/(e b £) « (56.1/e 6 ) /B[f] meV. 



(3.12) 



Cn agrees 2 ^ with the result of earlier work s 17 ' 18 obtained 
by a different method. 

This divergence in Cn derives from short- wavelength 
vacuum polarization and is present even for B = 0. To 
see this one may evaluate the Coulomb exchange correc- 
tion in free space (with a — > a/eb), using the instanta- 
neous photon and fermion propagators — 2ix a/(eb|k|) 
and iS(p) = (7^/(2 |p|), 



Y^v k iS(p + k)=<j-p£- log(CA 2 /p 2 ), (3.13) 



with momentum cutoff |k| < A and some constant C. 

This divergent correction causes infinite renormaliza- 
tion^ of velocity vq in the electron kinetic term oc vo^di; 
Eq. (|3.13l) agrees with an earlier result of Ref. [3(3. It van- 
ishes at p = but, for B ^ 0, turns into a nonvanishing 
energy gap, with p 2 — > 2eB = 2/£ 2 . Actually, this di- 
verging piece precisely agrees with that in Eq. (|3 . 1 1 1) . if 
one simply chooses the "Fermi momentum" A so that 
the Dirac sea accommodates the same number of elec- 
trons as in the B ^ case, N sca = A 2 /47r w N l /2tt£ 2 , 
i.e., A 2 « 2N L /£ 2 . 

Since such an infinite correction is already present for 
B = 0, it does not make sense to discuss the magnitude 
of the cutoff-dependent number Cat in Eq. (I3.11[) . The 
legitimate procedure is to renormalize vq by rescaling 



v = Z v Vn 



(3.14) 



, with 



VQ en regarded 



and put reference to the cutoff into Z, 
as an observable quantity. 

The renormalized velocity vj™ is defined by referring 
to a specific resonance. Let us take the — ► 1 resonance 
and choose to absorb the entire 0(V C ) correction at some 
reference scale (e.g., at magnetic field -Bo) into Z V7 i.e., 
we write 



by setting 



Ae 1 ^ = V2v^\ B /e 



(3.15) 



V2v 



■C N = 1 



8v e b 



A 2 

log -j, (3.16) 



where k 2 — (const.) eB$. The renormalized velocity then 
depends on B, or runs with B, 



ren i rcn 

v \b = v 



Bo - — \og(B/B ) 



(3.17) 



decreasing slightly for B > Bq; actually the correction is 
rather small (about 3% for B/Bq ~ 2 and e b ~ 5, with 
VQ Cn ~ c/300). With such B dependence in mind, we 
denote Vq cii \b as ^Jj™ and w£ en |s as w™ n from now on. 

The divergences in e£ xc = ej, — e n + Aek for all other 
resonances, as illustrated in Fig. 1 (a), are taken care of 
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FIG. 1: (a) Cyclotron resonance; circularly-polarized light 
can distinguish between two classes of transitions indicated 
by different types of arrows, (b) Momentum profiles of the 
many-body corrections Ae^J*""/ (\/n + 1 — y/n) in units of 
14 for n = 0,1,2 and 3. 



by this velocity renormalization. The finite corrections 
oc V c after renormalization then make sense as genuine 
observable corrections. In particular, for several intra- 
band channels L n — > L n+ i at filling factor v = 4nf + 2, 
direct calculations yield 

2<-l ^=6 
e k=o — 



2-l){^ cn - 0.264y c }, 
(V3-V2){o;^ n - 0.35814}, 
(V4- V3){^ on -0.419 Vc}, 
H-VZ){u™ -0.46414}, 



fc k=0 — 
e k=0 — 

elto U = 8 (VE - VZ) {u™ - 0.46414}, (3.18) 

where V c — ot/{e b l). The Coulomb corrections, shown 
numerically here, are analytically calculable. The excita- 
tion spectra ek in the hole band are essentially the same, 



-(rj+l)ii/=-(4ra+2) 
ln f =-(n+l) 



n+l<— n \v—4n+2 
k fif-Ti ' 



(3.19) 



reflecting the particle-hole symmetry. 

Figure 1 (b) shows some of the momentum profiles 



A 



-n+U 
k=0 



of Eq. (|3.9p , which, when integrated 



over £\p\, give Ae^+J^"/(Vn+ 1- y/n) in units of V c . It 
is clearly seen that the slowly decreasing high-momentum 
tails ~ (y/2/4)/(£ |p|) are responsible for the ultraviolet 
(UV) divergence and that the finite observable correc- 
tions oc V c are uniquely determined from the profiles in 
the low-momentum region I \p\ < 15. 

A look into the structure of the total current operator 
tells us that the optically-induced cyclotron resonance 
(for k = 0) in graphene is governed by the selection rule 
A|n| = ±1, in contrast to the "nonrelativistic" rule An — 
±1. In particular, there are two classes of transitions, 
(i) —n -> ±(n — 1) and (ii) ±(n — 1) — > n (with n > 1), 
which are distinguished^ by use of circularly-polarized 
light (oc A x ± iA y ); see Appendix B. 

As a result, graphene supports interband cyclotron res- 
onances. The lowest channels are open at v = — 2, with 



1 !/=-2 



c k=0 



__2 v=-1 



r k=0 



(V2 + l){< 
(V2 + 1)K 



0.12214}, 
0.15514}. (3.20) 
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FIG. 2: (a) Rescaled cyclotron-resonance energies as a func- 
tion of VB, with «5 cn w 1.13 x 10 6 m/s (at S = 10 T) and 
V c ~ 12 -y/ B [T] me V (et, ~ 5). The experimental data on 

eJ,.=o an d £k=o~ 1 /(v^2 + 1) ar e quoted from Ref. Q11, with er- 
ror bars inferred from the symbol size in the original data, 
(b) The same data plotted in units of ui c — \/2vo/t (with 
vo — > fo cn |-B=ioT) as a function of B. The dotted curve repre- 
sents a possible profile of the running of Vq c11 \b, normalized to 
1 at B = 10T, with ^ cn | s=10T « 1.13 x 10 6 m/s and e b « 5. 



Some other interband channels yield 

tlto 2 " = = 6 (V2 + 1){^ C " + 0.084 V c }, 
4to" 2 (V2 + V3){^ cn + 0.058 14}, 

4t<T 3 ^ 6 (V2 + V3)K C » + 0.H4 y c }, 

4^- 3 "== 10 (V2 + V3) {u;l cn + 0.04414}. (3.21) 

It is now clear that cyclotron resonance is best an- 
alyzed by plotting the rescaled energies e£.=o l\ s b\/W\ — 
s a -y/|a[| as a function of y/B or i3. The Coulombic many- 
body effect will be seen as a variation in the characteristic 
velocity «Q Cn [l + O(V^)] from one resonance to another, 
and a deviation of ui* c cn from the \f~B~ behavior would in- 
dicate the running of Wq™ with B. 

Figure 2 (a) shows such plots for some intra- and inter- 
band channels, using i>o ci1 |.b=iot ~ 1.13 x 10 6 m/s which 
fits the — > 1 resonance data, and, as a typical value, 
V c = a/(e b £) « 12 /B[f] meV (e 6 - 5). 

Actually, experiment^ has already observed a small 
deviation of the 1 : (1 + v2) ratio of £^=0 to c^Zq 1 well 
outside of the experimental errors under high magnetic 
fields B = (6 ~ 18) T; the data are apparently electron- 
hole symmetric, e k tTo ~ e k=o~ • Figure 2 (a) includes 



such data reproduced from Ref. [lj^. A small increase of 
Ug en in ej^TcT 1 / (V2+1), relative to e^o > is roughly consis- 
tent with Eq. (EES!]) which suggests a 0.122 V c /cj™ - 4% 
increase in Vg en (since V c /w™ n ~ 0.3). 

This feature is clearer from Fig. 2 (b), which plots the 
e ktTo ari d ^tlo /(v^ + data as a function of B in 
units of lu c = \/~2vq/£ cx \[B (with vq — > «o on |-B=ioT)- 
The deviation of the (—1 —> 2) resonance data is more 
pronounced. In the figure a dotted curve represents a 
possible profile of the running of Vq cu with B, and, es- 
pecially, the (—1 — > 2) data (with smaller error bars) 
suggests such running. 

It is too early to draw any definite conclusion from the 
present data alone, but the data is certainly consistent 
(in sign and magnitude) with the present estimate of the 
many-body effect. In this connection, let us note that 
an earlier experiment on thin epitaxial graphite 32 also 
observed the (0 — > 1) and (—1 — > 2) resonances, with 
apparently no deviation from the 1 : (1 + V2) ratio. This 
measurement was done under relatively weak magnetic 
fields B = (0.4 ~ 4) T, and it could be that a small devi- 
ation, under larger error bars, simply escaped detection, 
apart from the potential difference between thin graphite 
and graphene. 

More precise measurements of cyclotron resonance, es- 
pecially in the high B domain where the Coulomb inter- 
action becomes sizable, would be required to pin down 
the many-body effect in graphene. In this respect, the 
comparison between interband and intraband resonances 
from the same initial state, e.g., — n —> ±(n — 1) at 
v = 2 — An with n = 2, 3, would provide a clearer sig- 
nal for the many-body effect, with the influence of other 
possible sources reduced to a minimum. From Eqs. (|3.18p 
- (|3.21l) one can read off the variations in Vq 011 , 

Afl(-2->±l) ^ 6 0MV c /uj r c cn ~ 10%, 
AR(-3 -> ±2) " = « 10 0A0V c /uj r c cn ~ 12%, (3.22) 

which imply that a comparison of the (—2 — > ±1) res- 
onances and that of the (—3 — > ±2) resonances would 
find variations in vq, about 3 times larger than the ~ 4 % 
variation for ^=0" /(V^ + 1) vs e^Q 1 at v = —2. 

IV. CYCLOTRON RESONANCES IN BILAYER 
GRAPHENE 

In this section we consider cyclotron resonance in bi- 
layer graphene. In bilayer graphene the electrons are 
described by four-component spinor fields on the four in- 
equivalent sites (A,B) and (A',B') in the bottom and 
top layers, arranged in Bernal A'B stacking. Interlayer 
coupling 33 7i = ja'B ~ (0.3 — 0.4) eV modifies the in- 
tralayer linear spectra ±vq |p| to yield, in the low-energy 
branches |e| < 71, quasiparticles with a parabolic disper- 
sion^ They, in a magnetic field, lead to a particle-hole 
symmetric tower of Landau levels {L n } (n — 0±, ±1, ...) 
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with spectrum, 

e„ = s„w c ?7„((7i/w c ) 2 ), (4.1) 

1/2 

where T] n (x) = {(a n + x — \/x 2 + 2 a„ x + 1) /2} with 
a n = 2|n| — 1; see Appendix A. The sequence of low- lying 
levels is made clearer in the form 

e n = Sn uj hi ^\n\(\n\-l)/Uw), (4-2) 
with the characteristic cyclotron energy 

lo u = W 2 / 7l = 2v 2 Q /( ll f) ~ 5 B[T] meV, (4.3) 

where £ n (w) = {(l+a n w + \/l + 2 a n w + w 2 )/2} 1 2 and 

u> = (w c / 7 i) 2 = w c bi / 7 i - 0.01 SfT] < 1; £„(0) = 1. The 
high-energy branches |e| > 71 of the spectra give rise to 
another tower of Landau levels, with spectrum 

e+ = s„7i^nH, (4.4) 

where n = ±1, ±2, ■ • •. 

Note that £,„(w) = 1 + 0(e£>/ 7l ). As a result, e„ 
rises linearly with B at low energies |e n | < 71 and turns 
into a \f~B rise for |e„| 3> 71. Both e„ and e+ approach 
i-^/jnfwc for 1 71 1 — > 00, since the bilayer turns into two 
isolated layers at short wavelengths. 

In the bilayer there arise four zero-energy levels e„ = 
with n — (0±,±1) per spin. At one valley (say, K) they 
are electron levels with n = (0 + , 1) and, at another val- 
ley, they are hole levels with n = (0_,— 1); this feature 
is made explicit with a weak layer asymmetry, such as 
an interlayer voltage which opens up a (tunable) band 
gapj^i~— With a nonzero band gap, the zero-energy lev- 
els evolve into two quartets of nearly-degenerate levels 
(separated by the gap), i.e., "pseudo" -zero-mode levels, 
which are expected to support pseudospin wave o 34 ' 35 as 
characteristic collective excitations. 

For simplicity, we here turn off such a layer asymme- 
try as well as Zeeman splitting and the effect of trigo- 
nal warping (coming from 73 = 7 ab' < 71). In view of 
the small layer separation, we do not distinguish between 
the intralayer and interlayer Coulomb interactions. Each 
Landau level L n is thus treated as fourfold degenerate, 
except for the zero-mode levels (L 0+ ,Li) or (L _,L_i) 
which are fourfold degenerate at each valley. 

The effective Hamiltonian for the electrons in bilayer 
graphene takes a 4 x 4 matrix form which, for studying 
the properties of the low-lying levels, may be reduced to 
an approximate 2x2 form^ Actually, the bending of the 
spectrum e n with B is appreciable in the high-_B domain, 
B = (10 ~ 20) T, where cyclotron resonance in bilayer 
graphene has been studied experimentally. Accordingly 
we employ the full 4-component spinor description of the 
bilayer system; see Appendix A for details. 

The charge density p_ p (for each spin and valley) takes 
the same form as Eq. ()2.6[) . with g™ 1 replaced by 

= D k D n [fl k \M+p^^) f m-2M-2 
+(/3i 3) /3^+4 4) ^)/ifl- 1 '' n l- 1 ]; (4.5) 




FIG. 3: (a) Cyclotron resonance in bilayer graphene. 
(b) Momentum profiles of the bilayer many-body corrections 
AeS^7())„ + i - Tin) for n = 1,2, and 3, with v = 1.15 
and 71 = 3.5; real curves at B =10T and dashed curves at 
B =16T. A dotted curve refers to the profile of the monolayer 
(0 — > 1) resonance, (c) Resonance energies as a function of 
B; real curves, with Vo = 1.15, 71 = 3.5 and V c = 0; dot- 
ted curves, with -Co = 1.15, 71 = 3.8 and V c ~ 5.6 y/ B[T] 
meV (or e;, w 10). The experimental data with error bars are 
reproduced from Ref. l2ll . Note that the low-lying n = 2 spec- 
trum £2 (= the v = 4 curve) significantly deviates from the 
approximate spectrum e2 ~ splJ^. 1 with £2(1^) — > 1 (dashed 
line oc B) for B > 10T. 



see Appendix A for the coefficients } and D n . The 
sets (e n ,<7p ) at the K and K 1 valleys are related as 

en\K> = -e-n\K, 9 k p n \K> = Q^ltc (4.6) 

Actually, for zero band gap, are essentially the same 
at the two valleys since one further finds that 

ff£V =9?\k, 9%- 1 \k> = 9^\k,9^°-\k' = 9%° + \k, 

(4.7) 

for > 2 and \n\ > 2. 

One can now use the SMA formula (|3.9p to calculate 
the interlevel excitation energies e£ xc = — e a + Ae^~ a . 
The result applies to both valleys if one specifies the zero- 
mode levels accordingly. It is important to remember 
that for bilayer graphene the sum J2 n over filled levels 
involves two branches e_„ and in the valence band. 

Cyclotron resonance in bilayer graphene again obeys 
the selection rule 31 A|?i| = ±1; see Appendix B and 
Fig. 3 (a). The Coulombic corrections AeJ^Zg are diago- 
nal in spin and valley (while mixing arises for k ^ 0). The 
vacuum polarization effect again makes Ae£tg cutoff- 
dependent. 
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For renormalization let us first look into the B = 
case. One can construct the electron propagator and, as 
in the monolayer case, calculate the Coulombic quantum 
corrections. It turns out that not only vq but also 71 
undergo infinite renormalization and, rather unexpect- 
edly, the divergent terms are the same for both of them 
to 0(V C ) at least; they also coincide with the divergent 
term in the monolayer case; see Appendix C for details. 
To be precise, the divergences are removed, to 0{V C ) of 
our present interest, by rescaling 

v = Zv I ™, n = zir, (4.8) 

with a common factor Z. 

This scaling tells us how to carry out renormalization 
in the presence of a magnetic field B. Let us write, as in 
Eq. (|3.I8I) of the monolayer case, the excitation energy 
for the L n — > Lk transition in the form 

£k=o" = (skVk - s n vn) (V2v /i + c kn V c ), (4.9) 

with r] n = rj n (l/w). Note first that uo/71 = VQ en /j{ en 
is invariant under renormalization; it is therefore finite 
and does not run with B. Similarly, w = (uj c /-fi) 2 oc 
i v o en 1 7i en ) 2 B is invariant, and is linear in B. This means 
that t] n (l/w) remain unrenormalized and finite. Equa- 
tion (|4.9p then reveals a remarkable structure of the 
Coulombic corrections c kn : The divergent pieces are com- 
mon to all c kn and are removed by a single counterterm 
cx (Z- l)n5 CB 

Figure 3 (b) depicts the momentum profiles 7 2 [• • •] 
of Ae^^ 1 /{skTjk — s n i] n ) for some typical resonances. 
For comparison the profile for the monolayer resonance 
( Ae kio) mo °° is also included there. The gradually de- 
creasing high momentum tails, common to all, numer- 
ically demonstrate the validity of the scaling (|4.8p and 
Eq. (|4.9p . This further verifies that the leading loga- 
rithmic velocity renormalization is formally the same for 
both monolayer and bilayer graphene. 

For renormalization let us refer to a specific resonance, 
e.g., the — 3 — > — 2 resonance at v = —8, and define Vy en 
so as to absorb its entire 0(V C ) correction, 

wr = V2v™ n \ B /l=V2v /£ + c- 2 '- 3 V c . (4.10) 

One then has, for general n — > k channels, 

e k k to = (skVk ~ s nVn ) K on + Ac kn V c ). (4.11) 

Here Ac kn = c kn — c~ 2, ~ 3 are now free from the UV 
divergence and are uniquely fixed as genuine quantum 
corrections. In terms of the bilayer cyclotron frequency 
(cj bi ) rcn ee (u;™ n ) 2 /7j on , this also reads 

4to = (skCk - s n Cn) {(u hi Y cn + Ac fel Vw~V c }, (4.12) 

with Cn = VM(M ~ l)/€nH and w = (w bl ) ren /7f n - 

The quantum corrections Ac kn , unlike those of the 
monolayer case, are not pure numbers and, actually, are 



functions of y/w = w£ en /7™ n . This is seen if one notes 
that g kn are functions of y/w and £p so that c kn are func- 
tions of y/w and the cutoff Nl oc A 2 /(eB); the cutoff- 
independent corrections Ac kn thus depend on w alone. 
Let us set -/\ cn = 71 X 100 meV and vtf" = v a x 10 6 
m/s so that I /y/w = 7 ![ cn /a;™ n 2.75 G with G = 
7i/(7} 0v /B[Tj); G = 1 for 7l = 3.5 and v Q w 1.107 at 
B = 10T. It turns out that Ac fc '™, when plotted in G, 
behave almost linearly around G = 1. 

The way ^ en runs with B is determined from 

< n \ B = v r ™\ Bo +Sc- 2 <- 3 a/(V2e b ), (4.13) 

where 5c~ 2, ~ 3 = c~ 2 '~ 3 \g — c~ 2, ~ 3 \b - Numerically 
Sc~ 2, ~ 3 is nearly twice as large as the monolayer expres- 
sion -(V2/8) \og(B/B ) over the range 0.5 < B/B < 2 
around G = 1. The decrease in Vg en \B with B is larger 
in bilayer graphene and may amount to about 7% for 
B I Bo ~ 2 (and e& ~ 5). In this way, the renormalized 
velocity i>g cn is in general different, in magnitude and 
running with B, for monolayer and bilayer graphene; it 
reflects their low-energy features as well. 

We are now ready to look into some typical channels 
of cyclotron resonance. We use Eq. (I4.11[) and evaluate 
Ac 4 '" ee Ac numerically; for the bilayer the filling 
factor v = 4 (nf + 1) for m < —2 while v = 4nf for 
m > 1. For intraband channels one finds 

Ac ±i,-2 >^-4 Q.7270 + 0.5484 SG, 

Ac- 2 - 3 'V 8 0, 

Ac" 3 '" 4 v= = 12 -0.1521 - 0.0453 SG, 

Ac^'- 5 v= = 16 -0.2496- 0.0797 SG, (4.14) 

where SG = G - 1 with G = ^/{vq^/B]^]). Similarly, 
for interband resonances one obtains 

Ac 3 ^" 2 " = = i 0.3922 - 0.0023 SG, 

Ac 2 ^-' 3 " = = i 0.4794 + 0.0706 SG, 

Ac 2^- 3 ^v 8 0.3872 + 0.0552 SG, 

Ac 3^-4 ^=-12 Q.2961 + 0.00145 SG; (4.15) 

also Ac 4 ^- 3 "=~ 4 0.41 + ■ • • and Ac 4 ^ 3 " =~ 8 0.29 + - • •. 
These linearized expressions are numerically precise with 
errors of less than 3% over the range 0.3 < G < 1.5. The 
many-body effect is thus expected to be sizable in bilayer 
graphene. An effective variation in Ug on would amount to 
about 0.7V c /u} r c cn ~ 20% for e k tr o ~%=-4, and about -5% 

f° r e k=cT~ 4 |^=-i2> in comparison with e^ 2 _ ~ 3 | I/= _ 8 - 

As for experiment, Henriksen et al^ measured, via IR 
spectroscopy, cyclotron resonance in bilayer graphene in 
magnetic fields up to 18T. They observed intraband tran- 
sitions, which are identified with {ek=olf=4' e k=o\^=Si 
£^0^=12, e k=oU=i6} and the corresponding hole reso- 
nances listed in Eq. (|4.14p . together with an appreciable 
asymmetry between the electron and hole data. 



8 




FIG. 4: Experimental data of Ref . l2ll. reorganized in the form 
e k=o /( s kVk — s n T) n ) and plotted in units of lo c = \^2vo/£ 2 
(with do = 1-15 x 10 6 m/s). (a) Electron data, analyzed with 
«o = 1.15 and 71 = 3.5; for clarity the data points, originally 
at B=(W, 12, 14, 16) T, are slightly shifted in B. (b) The- 
oretical expectation according to Eq. (|4.11|l . with vq — 1.14, 
71 = 3.5 and V c /cj c w 0.24 (or e;, w 5.6). (c) Electron data, re- 
analyzed with 71 = 4. (d) Hole data, analyzed with vo = 1.02 
and 71 = 3.5. 



Figure 3 (c) reproduces the electron data of Ref. [2l|. 
There the real curves represent the resonance ener- 
gies (|4.11l) for V c = 0, with v ps 1.15 deduced from the 
v = 4 data and 71 taken to be 3.5, as supposed in Ref. I2H 
They poorly fit the v = 8, 12, 16 data. Unfortunately, in- 
clusion of the 0(V C ) corrections scarcely improves the fit, 
as seen from the dotted curves. 

The situation becomes clearer if one, in view of 
Eq. (|4.1ip . reorganizes the experimental data in the 
form £^=0 I ( s kVk — s n rj n ) and plots them in units of 



uj c = V2v /l (with v = 1.15 x 10 6 m/s). Figure 4 (a) 
shows such a plot for the electron data; for clarity the 
data points for different channels, originally at B=(10, 
12, 14, 16) T, are slightly shifted in B. It is to be con- 
trasted with Fig. 4 (b), which illustrates how each res- 
onance would behave with B, according to Eq. (|4.11l) . 
for V c ~ 10^/-B[T] meV (or w 5.6); in particular, the 
v = 8 curve represents the running of Vo° n |-B according 
to Eq. (|4.13[) . In Fig. 4 (a) the v = 8, 12, 16 resonances 
are apparently ordered in a way opposite to Fig. 4 (b), 
and an appreciable gap between the v = 4 resonance and 
the rest is not very clear. The v = 8, 12, 16 data show a 
general trend to decrease with B, consistent with possi- 
ble running of v r Q cn \B but at a rate faster than expected. 
It is rather difficult to interpret these features, but they, 
in part, could be attributed to possible quantum screen- 
ing^ of the Coulomb interaction in bilayer graphene such 
that Cb is effectively larger— for lower B. Note, in this 
connection, Fig. 4 (c) which shows that the same data 
may suggest a Coulombic gap for a choice 71 ~ 4 favored 
in Ref. l33l . We further remark that, in spite of an asym- 
metry in electron and hole data, the hole data shares 
essentially the same features; see Fig. 4 (d). 

No data are available for interband cyclotron resonance 
in bilayer graphene at present. They are highly desired 
because the comparison of interband and intraband res- 
onances from the same initial states would provide a 
clearer signal for the many-body effect. We record the 
ratios 

AR(-3 -> ±2) ^ 8 0.39 V c /ujl cn ~ 15%, 
AR(-4 -> ±3) " = w 12 0A5V c /u r c cn ~ 18%, (4.16) 

which imply that a close look into the (—3 — > ±2) res- 
onances and the (—4 — > ±3) resonances would find a 
sizable variation ~ 15% in Vy en . 

V. SUMMARY AND DISCUSSION 

Graphene supports charge carriers that behave as 
Dirac fermions, which, in a magnetic field, lead to a char- 
acteristic particle-hole symmetric pattern of Landau lev- 
els. Accordingly, unlike standard QH systems, there is a 
rich variety of cyclotron resonance, both intraband and 
interband resonances of various energies, in graphene. 

In this paper we have studied many-body corrections 
to cyclotron resonance in graphene. We have constructed 
an effective theory using the SMA and noted that genuine 
nonzero many-body corrections (not due to fine splitting 
in spin or valley) derive from the quantum fluctuations 
of the vacuum (the Dirac sea). Such quantum correc- 
tions arc intrinsically ultraviolet divergent and, as we 
have emphasized, it is necessary to carry out renormaliza- 
tion of velocity vq (and, for bilayer graphene, interlayer 
coupling 71 as well) to determine the many-body cor- 
rections uniquely in terms of physical quantities. As a 
result, the observable intralayer and interlayer coupling 
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strengths Wq 6 " oc 7jj en and 7J 011 in general run with the 
magnetic field B. 

Experimental data on cyclotron resonance generally 
have sizable error bars, which make a clear identification 
of the many-body effect difficult. In this respect, we have 
presented a way to analyze the data, as in Fig. 2 (b) and 
Fig. 4, with the effect of renormalization properly taken 
into account. 

For monolayer graphene a piece of data^ which com- 
pares some leading interband and intraband resonances 
is apparently consistent with the presence of many-body 
corrections roughly in magnitude and sign, and also in 
the running of vfj en with B. 

For bilayer graphene the existing data are only for in- 
traband resonances and are rather puzzling, as discussed 
in Sec. IV. They generally appear to defy good fit by 
theory but certainly suggest nontrivial features of many- 
body corrections, such as running with B. 

More precise measurements of cyclotron resonances are 
highly desired. Of particular interest are experiments 
which compare interband and intraband resonances from 
the same initial states, as listed in Eqs. (I3.22[) and (|4.16l) . 
which would clarify the many-body effect with minimal 
uncertainties. 



The corresponding eigenfunctions for n 
take the form 



±2, ±3, 



(\ 



/3i 2) |M-2) 

/?i 3) IM-i> 



\ 



V/?i 4) IH 



(A3) 



r'n 1 r'n 



1 \n 

~7 



,D n = -7= 



'|n|(|n|-l- e ' 2 ) 



V2 I/ M(M-l)-e' 4 



(A4) 



where only the orbital eigenmodes are shown using the 
standard harmonic-oscillator basis {!«)}• These expres- 
sions for \& n are equally valid for both the low- and high- 
energy branches e n and e+ of Landau levels, depending 
on e' one employs. 

The zero-energy eigenmodes are given by 



* 0+ = (|0), 0,0,0)*, 
* 1 = J D 1 (|l),0,-(l/ 7 ')|0),0) t 



(A5) 
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Appendix A: Landau levels in bilayer graphene 

This appendix summarizes the effective Hamiltonian 
and its eigenfunctions for bilayer graphene in a magnetic 
field B. The bilayer Hamiltonian with interlayer coupling 
71 = "fA'B is written, at one (K) valley, as^ 



H 



hi 



( v nt^ 
u n 
« nt 
\ i> n 71 j 



71 




(Al) 



and 11+ = U-a + illy, with A t -t B(-y,0). 
The energy eigenvalues obey the equation 



(|n| 



1 



e' 2 )(\n\ 



e )-7 (O 



0, 



(A2) 



where e' = e n /ui c , 7' = 7i/u; c and cj c = ^/2vo/£- This 
leads to the two branches of spectra (e„, e+) in Eqs. (|4.1[) 
and (I4.4p . In particular, zero energy e„ = is possi- 
ble for |n| = or |n| = 1 while \e±i\ > 71. A weak 
interlayer voltage \{5V) diag[l, —1, —1, 1], added to H h \ 
reveals that the zero modes actually have n = 0+ and 
n = 1 for SV > 0. 



with Lq = (1 + 1/7' 2 )- 1 / 2 . 

At another (K') valley the Hamiltonian is given by 
Eq. (|Alj) with vq — > — vq and acts on a field of the form 
^ K' = (ipB 1 , ipA, ipB, tpA'Y- Accordingly one finds that 

e n \K' = -e-n\K,D n \ K * = D- n \ K ,f3^\ K , = /3^1\k, 
A 3) \k' = -/3^k,/?i 4) k'=-^k- (A6) 
The zero-energy levels now have n — 0_ and n = — 1. 

Appendix B: Coupling to current 

Consider a weak time-varying vector potential 
(A x (t), A y (t)) coupled to the total current in graphene. 
For the effective Lagrangian L$ in Eq. (|3.4[) this yields 
coupling of Ai to = / d 2 x $ kn of the form 



which acts on an electron field of the form ^> k = ^ A ~ 
(ipA, V'B', ipA', iPbY in obvious notation; n = — iIL y 



-i^fp- Q d n {A^ Q 



V2 
+h.c, 



±(n-l),-n 



n,±(n-l) 



k=0 



} 

(Bl) 



where d n = ±1/2 for n > 2 and d 1 = il/%/2; A = A x - 
iA y ] po = l/(2ir£ 2 ). The cyclotron resonance thus obeys 
the selection rule A|n| = ±1. In particular, the —n^r 
±(n — 1) transitions and the ±(n — 1) — » n transitions 
(n > 1) are distinguished^ by use of circularly-polarized 
light oc A x ± iA y . 

Equation (IJ3L|) (with n > 2) applies to the case of bi- 
layer graphene as well if one sets oj c — > uj^ 1 , A^ —> —A', 
d n = \/n — 1 for n > 3 and d 2 — \/2, apart from terms 
ofO((^/ 7 i) 2 ). 
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Appendix C: propagators 

In this appendix we derive the electron propagator for 
bilayer graphene in free space. Let us set II — > p x — ip y 
in H hl of Eq. (|Aip and consider the propagator 

(*(a:)*V)) = (x\l/{id t - H hi )\x') (CI) 

with |x) = |x) \t). We divide the 4x4 matrix H hl into a 
2x2 block form and invert (idt — H hl ). In Fourier (p, w) 
space the propagator reads, in 2 x 2 block form, 

( W t) n = A| w ( w 2_ v 2 p 2_ 7 2 )+7it; 2 p(TiP | 5 

(*$t) 12 = l-voP^ 2 -v$p 2 + W7i<Ti), 

= » ~ - «§ P 2 + w 7i cti), (C2) 

where D = (« 2 p 2 - w 2 ) 2 - u; 2 7l 2 = (a; 2 - E 2 + )(lu 2 - E 2 _) 
with E± = ^7^/4 p 2 ± 7i/2; P = cr + + pcr_ and 



Pcri-P = (p t ) 2 cr + +p 2 a-_ withp = p x ~ip y , =p x + ip y 
and <r± = (<Ti ± io2)j2. 

This leads to the instantaneous propagator 
/(dw/27r)(**t) j 

/**t\| _ 1 /^71^1^/P 2 2 V0 P\ fr „s 



with V p = V7?/4 + « 2 P 2 . 

To calculate the Coulomb exchange correction one may 
replace, in Eq. fl3.13p . iS(p) by this propagator. Note 
that vq P/(2T> P ) approaches, for p — > oo, the mono- 
layer propagator iS(p) = <JiPi/(2 |p|) (apart from an 
inessential mismatch ai — > — (72 in notation). As a re- 
sult, setting iS(p) —> vo P/ (2V p ) for vq and iS(p) — >• 
7i/(42? p ) for 71 and carrying out the k integration, as in 
Eq. (|3.13p . yield the same amount of logarithmic diver- 
gence ~ (a/8efc) log A 2 as in the monolayer case; it thus 
renormalizes v and "fi simultaneously as in Eq. (14.81) . 
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